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We analyze the effect of blood flow through large arteries of peripheral circulation on heat transfer 
in living tissue. Blood flow in such arteries gives rise to fast heat propagation over large scales, which 
is described in terms of heat superdiffusion. The corresponding bioheat heat equation is derived. 
In particular, we show that under local strong heating of a small tissue domain the temperature 
distribution inside the surrounding tissue is affected substantially by heat superdiffusion. 



I. INTRODUCTION. LIVING TISSUE AS A 
HETEROGENEOUS MEDIUM 

Blood flowing through vessels forms paths of fast heat 
transport in living tissue and under typical conditions 
it is blood flow that governs heat propagation on scales 
about or greater than one centimeter (for an introduction 
to this problem see, e.g., ||)- Blood vessels make up a 
complex network being practically a fractal. The larger 
is a vessel, the faster is the blood motion in it and, so, 
the stronger is the effect of blood flow in the given vessel 
on heat transfer. Blood flow in capillaries practically 
does not affect heat propagation whereas blood inside 
large vessels moves so fast that its heat interaction with 
the surrounding cellular tissue is ignorable 0. Thus, 
there should be vessels of a certain length that are 
the smallest ones among the vessels wherein blood flow 
affects heat transfer remarkably. The value of can be 
estimated as |3| (see also I, |): 
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where D = k/ {cp) is the temperature diffusivity of the 
cellular tissue determined by its thermal conductivity k, 
specific heat c, and density p, the value j is the blood 
perfusion rate (the volume of blood going through tis- 
sue region of unit volume per unit time), and the fac- 
tor Ln, ~ \n.{l/a) is logarithm of the mean ratio of the 
individual length to radius of blood vessels forming pe- 
ripheral circulation. For the vascular networks made up 
of the paired artery and vein trees where all the vessels 
are grouped into the pairs of the closely-spaced arter- 
ies and veins with opposite blood currents the coefficient 
/ ^ Ln^^^ accounts for the counter-current effect [ p^ . 
For the vascular networks where the artery and vein trees 
are arranged independently of each other the factor / 
should be set equal to unity, / = 1. In particular, for the 
typical values of the ratio l/a ^ 20-40 [Q, the thermal 
conductivity k ~ 7 • 10^'^ W/cm-K, the heat capacity 
c ~ 3.5 J/g- K, and the density p ^ 1 g/crn^ of the tissue 
as well as setting the blood perfusion rate j ^ 0.3min~^ 
from (^ we get the estimates ~ 4 mm and Ln ~ 3 — 4 . 
In the mean-field approximation the effect of blood 




vein node 



FIG. 1: Schematic illustration of the effect of blood flow 
through the vein tree on heat diffusion imitated by random 
walks. The figure shows trapping of a random walker because 
of getting the internal points of a large vein after passing the 
vein node. 



flow on heat transfer is reduced to the renormalization 
of the temperature diffusivity, D — > -Dcff, Jp| and the 
appearance of the effective heat sink fj [0, |^|^ in the 
bioheat equation: 

^ = V(DeffVr)-/j(T-rj + gT. (2) 

Here T(r, t) is the tissue temperature field averaged over 
scales about the parameter Ta is the blood temper- 
ature insider the systemic circulation arteries, and the 
summand qT{v,t) called below the temperature gener- 
ation rate is specified by the heat generation rate q as 
It = 9/(cp). The renormalization of the temperature 
diffusivity is mainly determined by the blood vessels of 
lengths about and due to the fractal structure of vascu- 
lar networks the renormalization coefficient F — D^s/ D 
is practically a constant of unity order, > 1 0. 

Let us imitate the temperature evolution in terms of 
random walks whose concentration is {T — Ta). Then the 
part of the vein tree made up of vessels whose lengths 
exceed or are about the scale forms the system of traps. 
In fact, blood streams going through the vein tree merge 
into greater streams at the nodes (Fig. |l]). Therefore 
an effective random walker after reaching the boundary 
of one of these veins inevitably will be moved by blood 
flow into the internal points of large veins. Then, due to 
relatively fast blood motion inside these veins it will be 
carried away from the tissue region under consideration, 
which may be described in terms of the walker trapping 
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or, what is the same, the heat sink [^j. Since the mean 
distance between these veins is determined mainly by the 
shortest ones, i.e. by the veins of length the mean time 
during which a walker wanders inside the cellular tissue 
before being trapped is |^ 

-^-^L^- (3) 

In obtaining the given expression we have assumed the 
vascular network to be embedded uniformly in the cellu- 
lar tissue, so the tissue volume £^ falls per one vein (and 
artery, respectively) of length Whence it follows, in 
particular, that the rate at which the walkers are being 
trapped by these vein, i.e. the rate of their disappear- 
ance is estimated as 1/t, leading together with expres- 
sion (|l|) to the heat sink of intensity fj in the bioheat 
equation (||). The characteristic spatial scale of walker 
diffusion in the cellular tissue before being trapped is 
~ Vi>r, i-e. 

It ^ \ ^Icm. (4) 

V ]J 

The scale gives us also the mean penetration depth 
of heat penetration into the cellular tissue from a point 
source or, what is the same, the widening of the temper- 
ature distribution caused by heat diffusion in the cellular 
tissue. It is the result obtained within the mean-field 
approximation. 

Beyond the scope of the mean- field theory we meet sev- 
eral phenomena. One of them is the temperature nonuni- 
formities caused by the vessel discreteness Q which can 
be described assuming the heat sink in equation (|^) to 
contain a random component j^. Another is fast heat 
transport over scales exceeding substantially the length 
caused by the effect of blood flow through the artery 
tree. The latter phenomena is the main subject of the 
present paper. 

II. FAST HEAT TRANSPORT WITH BLOOD 
FLOW THROUGH LARGE ARTERY TREE 

Let at a certain time a random walker wandering in 
the cellular tissue get a boundary of a large artery, i.e. 
an artery of length exceeding t^. It should be noted 
that such an event is of low probability and cannot be 
considered within the standard mean-field approximation 
because the relative number of large arteries is small. 
Due to the direction of the blood motion from larger 
arteries to smaller ones as well as the high blood flow 
rate in the large arteries the walker will be transported 
fast to one of the arteries of length (Fig. ^) . The blood 
flow rate in small vessels is not high enough to affect the 
walker motion essentially and it has inevitably to leave 
this artery and wander in the cellular tissue until being 
trapped by the veins of length ^y. Thereby a certain 
not too large number of random walkers generated, for 




FIG. 2: Schematic illustration of the anomalously fast heat 
propagation caused by blood flow through large artery tree in 
terms of random walks. The figure shows the effective walker 
motion with blood from a large artery to a small one where 
the walker leaves it wandering in the cellular tissue. 



example, inside a cellular tissue neighborhood of a point 
r can be found during the time r inside a cellular tissue 
neighborhood of a point r' at a distance much larger than 
^t,, i.e. |r — r'l ^ The given effect may be regarded as 
anomalously fast heat diffusion in living tissue, i.e. heat 
superdiffusion. 

Dealing with heat transfer in living tissue we may con- 
fine our consideration to the peripheral vascular networks 
typically embedded uniformly into the cellular tissue, at 
least at the first approximation j|] . The latter statement 
means, in particular, the fact that for a fixed peripheral 
vascular network the vessel collection comprising all the 
arteries of length I meets the condition of the volume 

approximately falling per each one of these arteries. 
Therefore as is seen in Fig. ^ a walker going into a large 
artery of length I at initial time during the time r before 
being trapped by the veins can be found equiprobably at 
each point of the given artery neighborhood of size I. In 
other words, this walker makes a large jump of length I 
that exceeds substantially the mean-field diffusion length 

In what follows we will analyze the temperature distri- 
bution averaged over all the possible realizations of the 
vascular network embedding. This enables us to regard 
a walker entering a large artery of length I as a random 
event whose probability is independent of the walker ini- 
tial position. Thereby, the probability Pi for a walker to 
make a large jump over the distance I is also indepen- 
dent of the spatial coordinates r. It should be noted that 
for a fixed realization of the vascular network embedding 
the probability P/ depends essentially on the spatial co- 
ordinates r and the heat transfer in living tissue on large 
scales has to exhibit substantial dependence on the spe- 
cific position of the heat sources. 

Now we estimate the value of Pi assuming the heat 
sources to be localized inside a domain Qc of size L. Two 
different factors determine the value of Pi. First, it is the 
process of walker trapping by an artery of length I > C 
going through the domain Qc- If I < C blood flow in this 
artery has practically no effect on heat diffusion. On the 
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average a random walker during the time r travels the 
distance £t in the cellular tissue until being trapped by 
the veins of length £„. So for a walker to enter this artery 
and, thus, to leave the domain Qc with blood flow in it 
the walker, on one hand, should be located at initial time 
inside a cylindrical neighborhood Qi of the given artery 
whose radius is about ix and the volume is Ci"^. On the 
other hand, it has to avoid being trapped by the veins 
of length £y. The probability of the latter event is about 
(£„/£t)^- Indeed, a vein of length may be treated as a 
trap of cylindrical form. Thereby in qualitative analysis 
the walker trapping can be described in terms of two- 
dimensional random walks in the plane perpendicular to 
the artery under consideration where the trapping veins 
are represented by small circular regions. Their density 
is about l/if, which directly leads to the latter estimate. 
Therefore the total number of walkers leaving the domain 
Qc with blood flow through the given artery per unit 
time is 



■ al ■ [T - Ta) = ^C{T - Ta) . (5) 



In obtaining we have taken into account expres- 
sion . Since the trapped walkers spread uniformly over 
a region of size / the resulting density of the walker tran- 
sition rate to a point r spaced at a distance about I from 
the domain Qc is 



(6) 



It should be noted that the transition rate gi{r), as it 
must, does not depend on the local value of blood perfu- 
sion rate. 

At the second step we should average the obtained 
transition rate gi (r) over the possible realizations of the 
vascular network embedding. Let us adopt a simplified 
model for the vascular network shown in Fig. [^a where 
the vessel lengths In and In+i of the neighboring hierar- 
chy levels n and n + 1 are related as /„ — 2ln+i- Figure ||& 
demonstrates a more adequate model for the peripheral 
artery tree which, however, within the framework of the 
present qualitative analysis may be reduced to the for- 
mer one by combining three sequent two-fold nodes into 
one effective four-fold node at all the levels. In this case 
the cubic domain of volume 1^ falls per each artery of 
level n. Let us now consider individually three charac- 
teristic forms of the domain Qc, a ball or a cube of size 
£ (d = 3), an infinitely long cylinder of radius £ (d = 2), 
and a plane layer of thickness C (d = 1). For the ball or 
the cube, i.e. a region bounded in three dimensions the 
probability that an artery of level n passes through the 
domain is about 



P, 



(3e) 





FIG. 3: Models of the peripheral artery network embedding 
into the cellular tissue, (a) four-fold node model used in the 
present analysis and (6) a more realistic dichotomic artery 
tree uniformly embedded into a cellular tissue domain M. In 
the qualitative description of heat transfer both the models 
lead to the same result. 



For the infinitely long cylinder 

,(2e) / £ 



and for the plane layer Pi ^''' ^ 1. Multiplying gi{r) by 

the corresponding values of P^'^'' we get the result of av- 
eraging the walker transition rate gi (r) over the possible 
realizations of the vascular network embedding. The ob- 
tained result is written as 
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where the value d actually plays the role of the dimension 
of the space inside which the temperature field can be 
considered. At the next step we should sum the terms (0) 
over all the levels of the large artery tree. However due 
to the strong increase of the terms (0), (5;(r)) 2"(^+'*), 
as the level number n increases the arteries of length 
I ^ r mainly contribute to the value of {gi{r)). So the 
term describing the fast heat transport with blood flow 
through large arteries from the domain Qc (located near 
the origin of the coordinate system) can be written as 
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where Ai is the region containing the peripheral vascu- 
lar network as a whole and the integration in the three 
dimensional space over the domain Qc allows for all its 
three considered types. 

Expression (||) together with the mean-filed bioheat 
equation (|^) enables us to write the following equation 
governing the anomalous heat transfer in living tissue 

dt LnJM \y-y'Y+£1 

(9) 

where we have added directly the value ly in order to 
cut off the spatial scales smaller than the length and 
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ignored the difference between the effective temperature 
diffusivity and the true one of the cellular tissue. Equa- 
tion ^ is the desired governing equation of the anoma- 
lous fast heat diffusion in living tissue for the averaged 
temperature field. It should be noted that the second 
term on the right-hand side of equation ^ depends 
weakly on the blood perfusion rate. Therefore, for the 
nonuniform distribution of the blood perfusion rate j(r, t) 
equation holds also. 

Anomalous temperature distribution under local 
strong heating 

Hyperthermia treatment as well as thermotherapy of 
small tumors of size about or less 1 cm is related to local 
strong heating of living tissue up to temperatures about 
45 °C or higher values. In this case the tissue region 
heated directly, for example, by laser irradiation is also 
of a similar size. Due to the tissue response to such strong 
heating the blood perfusion rate can grow tenfold locally 
whereas in the neighboring regions the blood perfusion 
rate remains practically unchanged The feasibility 
of such nonuniform distribution of the blood perfusion 
rate may be explained applying to the cooperative mech- 
anism of self-regulation in hierarchically organized active 
media Therefore in the region affected directly 

the blood perfusion rate jmax can exceed the blood per- 
fusion rate jo in the surrounding tissue substantially. In 
this case the characteristic length of heat diffusion into 
the surrounding tissue is about 

y J Jmax 



giving us also the minimal size /Cmim of the region wherein 
the tissue temperature increase (Tmax — 7a) is mainly lo- 
calized. In the neighboring tissue the blood perfusion 
rate keeps up a sufficiently low value jo i which makes the 
heat propagation with blood flow through large arteries 
considerable. Indeed, let us estimate the temperature 
increase caused by this effect using the obtained equa- 
tion (|^). The temperature increase T(r) — at a point 
spaced at the distance r > C from the region (of size 
C) affected directly, i.e. the tail of the temperature dis- 
tribution is mainly determined by the anomalous heat 
diffusion and, so, is estimated by the expression 



rW-r.^|^^%^(T„_-T,) . (11) 



As seen from (|T^) for a sufficiently local and strong heat- 
ing of the tissue, i.e. when C ^ and jmax ^ jo the 
temperature increase at not too distant points such that 
r ^ C can be considerable. Otherwise the anomalous 
heat diffusion is ignorable. 
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